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Direct photon detection is experimentally implemented to measure the squeezing and purity of 
a single-mode squeezed vacuum state without an interferometric homodyne detection. Following a 
recent theoretical proposal [arXiv quant-ph/0311119 , the setup only requires a tunable beamsplitter 
and a single-photon detector to fully characterize the generated Gaussian states. The experimental 
implementation of this procedure is discussed and compared with other reference methods. 
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I. INTRODUCTION 

Squeezed states of light play an important role in the 
development of quantum information processing with 
continuous variables 0, where the information is en- 
coded in two conjugate quadratures of an optical field 
mode. These states may for example be used as a main 
resource for quantum cryptographic protocols (see 0, Q 
and references therein). They may also serve as an en- 
tanglement source since combining two squeezed states at 
a beamsplitter creates an entangled two-mode squeezed 
state such as those required for quantum teleportation 
or dense coding 0- In addition, squeezing has been 
shown to be an irreducible resource for realizing an arbi- 
trary linear canonical transformation 6]. 

Any attempt to process squeezed states in quantum 
communication or computation systems will necessarily 
face the problem of characterizing these states. A pos- 
sible complete description of a general quantum state 
is obtained by reconstructing its Wigner function using 
quantum tomographic procedures 7,8,9]. Alternatively, 
for a Gaussian state, which is fully described by its first 
and second order moments, another complete character- 
ization is provided by the mean values of the conjugate 
quadratures x and p together with the associated covari- 
ance matrix 7. From this, one may compute various rele- 
vant parameters such as the maximum observable squeez- 
ing 01 or the degree of purity 0, 0] . 

In this paper, we follow an idea originally due to two of 
the present authors that consists of measuring the squeez- 
ing and purity of a Gaussian state without homodyne de- 
tection, that is, without any strong local oscillator beam 
providing a phase reference 0. The suggested setup 
only relies on beamsplitters and single-photon detectors. 
It generally requires a joint measurement of two copies 
of the Gaussian state, but single-copy measurements suf- 
fice if it is a priori known that the mean values of the 
quadratures vanish. Thus, in the latter case, which ac- 
tually applies to all quantum information schemes based 
on squeezed vacuum states, no interferometric stability 
is required to determine the squeezing and purity, unlike 
with homodyne detection schemes. 

Hereafter, we will focus on the important case of 



a single-mode squeezed vacuum state and discuss the 
experimental feasibility and relevance of this photon- 
counting characterization procedure. Some useful no- 
tations to describe a squeezed vacuum beam are intro- 
duced in section [HI Section ITTT1 then presents the exper- 
imental setup together with two classical and homodyne 
measurement procedures that are used as a reference to 
characterize the generated squeezed vacuum states. In 
Section lTVl we briefly review the photon-counting charac- 
terization method applied to the special case of a single- 
mode Gaussian state (more details can be found in [l3|h 
In Section we present the experimental results of this 
characterization method, while Section fVTl discusses the 
constraints on the global detection efficiency that are put 
by this method. Numerical simulations are used to illus- 
trate the photon-counting method for values of the global 
detection efficiency that are presently unreachable in the 
experiments. 



II. SQUEEZED VACUUM DESCRIPTION 

Theoretically, a general Gaussian state with zero mean 
values of quadratures is fully characterized by its covari- 
ance matrix 7, which comprises the second moments of 
the conjugate quadratures x = a + a^ and p = (a — a')/i, 
with [x,p] — 2i. For states with zero mean values of 
quadratures, the covariance matrix 7 can be expressed 
as follows, 
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In order to determine the squeezing and purity, we only 
need to measure the two invariants of the covariance ma- 
trix, namely the trace Tr(7) and determinant det(7). 

From a more physical point of view, one can use the 
fact that the most general single-mode Gaussian state 
with (x) = (p) = can be expressed as a squeezing 
operator applied to a Gaussian thermal state |12L Il5| . 
Translated into an optical setup, this is implemented by 
two simple linear amplifiers as depicted in Fig. ^ : a 
phase-insensitive amplifier of gain H followed by a phase- 
sensitive amplifier of gains G and 1/G (in the following 
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FIG. 1: General single- mode Gaussian state generation. H 
denotes the intensity gain of the phase-independent amplifier 
and G stands for the intensity gain of the phase-dependent 
amplifier. 
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we take G, H > 1). In other words, the physics of the 
optical parametric amplifier (OPA) can be modeled by a 
"black box" squeezer which is parametrized by H and G. 
These two parameters are equivalent to the two phase- 
insensitive parameters Tr(7) and det(7) of the Gaussian 
state generated by the black box from the vacuum. 

Let us describe the transformation effected by the am- 
plifier depicted in Fig. QJ First, one can express the con- 
jugate quadrature variables at the output of the OPA 
as 



•'nul = — = ! \ U.l\. it , -'r \ II - U aJlc ), (2) 

Pout = VG(V~Hp vac - VH - lpanc), (3) 

where we chose x ou t (Pout) as the squeezed (anti- 
squeezed) quadrature, and x vac and x anc denote the vac- 
uum and ancilla quadratures at the input ports of the 
total amplifier, respectively. One can then express the 
variances of the squeezed and anti-squeezed quadratures 
at the output of the OPA as 



(2.ff- 1)/G, 
(2H- 1) G. 
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The trace and determinant of the covariance matrix read 
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This system of equations can be inverted, and the 
squeezed and anti-squeezed variances can be expressed 
in terms of the trace and the determinant of 7, 
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Finally, the purity V — Tr[p 2 ] of a mixed state p is, for 
any single-mode Gaussian state, directly linked to the 
average photon number of thermal noise n = H — 1 : 



V = 



1 



1 



(9) 



2n + 1 2H-1 
Equivalently, in terms of the covariance matrix, we have 

1 



V = 



Vdet( 7 ) 
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FIG. 2: Simplified experimental setup. SHG : second har- 
monic generation module; OPA : degenerate optical paramet- 
ric amplifier; APD : avalanche photodiode photon counting 
module. 



III. REFERENCE CLASSICAL AND 
HOMODYNE CHARACTERIZATION METHODS 

A new scheme for pulsed squeezed light generation has 
recently been developed Q and will be used here to 
compare the photon-counting characterization method to 
standard methods. The experimental setup is depicted in 
Fig. |21 The initial pulses are obtained from a titanium- 
sapphire laser delivering nearly Fourier-transform limited 
pulses centered on 846 nm, with a duration of 150 fs, a 
typical energy of 40 nJ, and a repetition rate of 780.4 
kHz. These pulses are frequency doubled in a single pass 
through a thin (100 /xm) crystal of potassium niobate 
(KNbOa), cut and temperature-tuned for non-critical 
type-I phase-matching. The second harmonic power is 
large enough to obtain a significant single-pass paramet- 
ric gain in a similar KNb03 crystal used in a type-I spa- 
tially degenerate configuration. 

The squeezed beam can then be directed onto two dif- 
ferent detection modules using a removable mirror: 

(i) Homodyne detection module : the squeezed vacuum 
beam interferes with the local oscillator beam (LO) in a 
balanced homodyne detection setup. A main feature of 
our experiment is that all the processing is done in the 
time domain, not in the frequency domain. For each 
incoming pulse, the fast acquisition board samples one 
value of the signal quadrature in phase with the local 
oscillator [l4|. 

(ii) Photon counting module : the squeezed vacuum 
beam is transmitted via a beam splitter of tunable trans- 
mittance T and then passes through a spatial filter (made 
of two Fourier-conjugated pinholes) and a 3 nm spectral 
filter centered at the laser wavelength, before being de- 
tected by a silicon avalanche photodiode (APD). 

To start the characterization procedure, a basic mea- 
surement is to monitor the classical amplification and 
de-amplification of a probe taken from the fundamen- 
tal beam. This is easily done by direct detection of a 
probe beam averaged power on a photodiode. Setting 
the relative phase between the probe and the second har- 
monic pump beam allows to tune the classical gain from 
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the minimum de-amplification intensity gain <7 m in to the 
maximum gain £/ max . The measurement of the classical 
gains Qiain and £/ max gives an estimate of G and H , 



G = \J Gmax/G n 
^ \/^max^min ■ 



(ii) 

(12) 



The experimental results of the squeezed vacuum char- 
acterization for different values of the pump power are 
shown in Figs. [3] and marked as "Classical" (black 
disks). 

Following the principle of quantum tomography, a 
powerful approach is to completely characterize the 
squeezed vacuum by conjugate quadratures homodyne 
measurements. The time-resolved balanced homodyne 
detection allows to measure the squeezed and anti- 
squeezed quadrature variances Vh om ,min and Vh oro ,max- 
Imperfections and losses in this detection are modeled 
by a beamsplitter of transmission ryi lom (in intensity). 
The procedure to measure the detection efficiency is well 
established from squeezing experiments [l6j . and it can 
be cross-checked by comparing the classical parametric 
gain and the measured degree of squeezing. We note 
the homodyne detection efficiency % m = VtVhVd — 
0.76 ± 0.01, where the overall transmission r\T = 0.92, 
the mode-matching visibility rj H = 0.935, and the detec- 
tors efficiency t]d = 0.945 are independently measured. 
Given this efficiency, one can correct for losses and deduce 
the squeezed and anti-squeezed quadrature variances at 
the output port of the OPA, namely 
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This allows the full characterization of the state param- 
eters (Tr(7), det(7)) or the OPA parameters (G, H) fol- 
lowing the above formulae. The experimental results of 
this second characterization method are also displayed 
in Figs.|21and21 marked as "Homodyne" (gray squares). 
As one can notice in Figs. |3]and01 for high pump pow- 
ers the "homodyne" and "classical" values do not well 
overlap within their respective error bars. A main rea- 
son for this is that the "black-box" model developed 
above is basically a single-mode model, and thus suffers 
from fundamental limitations, while for high pump pow- 
ers the physics involved in parametric de-amplification 
are known to fall into a multi-mode regime |17| . 

In the following, we will use these "classical" and "ho- 
modyne" methods as references to check the validity of 
the photon-counting characterization. 



IV. PHOTON-COUNTING 
CHARACTERIZATION METHOD 

Let us briefly introduce the methods implemented here 
for measuring the properties of a Gaussian state by pho- 
ton counting, which are derived from the original proce- 
dure presented in |13J . We will restrict our attention to a 
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FIG. 3: Parametric gains G, H and purity V versus average 
pump power at 423nm. The solid line corresponds to a fit on 
the classical results according to plane wave theory. "Classi- 
cal" stands for the classical probe gain measurements (black 
disks), "Homodyne" stands for the balanced homodyne de- 
tection variance measurements (gray squares). The bounds 
inferred from the photon counting method, to be described in 
section Hvl are indicated by the two dashed lines (see further 
explanations in section IVt . 



single-mode Gaussian state with zero coherent displace- 
ment, (x) = (p) = 0, which is the case in the present 
experiment depicted in Fig. The squeezed vacuum 
mode impinges on a beam splitter with tunable trans- 
mittance T before being measured by an avalanche pho- 
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FIG. 4: Trace Tr(7) and determinant det(7) versus average 
pump power at 423nm. The annotations are the same as in 
figED F° r ease °f viewing, the trace values obtained from 
the photon counting method (gray diamonds) are linked by a 
dash-dotted line in Fig. (a). In Fig.(b), the two dashed lines 
indicate the limits on the determinant knowledge obtained 
from the photon counting method (see further details in sec- 
tion El • 



todiode that is sensitive to single photons and can re- 
spond with two measurement outcomes, either a click or 
a no-click. This detector with overall detection efficiency 
?7apd can be modeled as a beam splitter with transmit- 
tance ?7apd followed by an ideal detector that performs 
a dichotomic measurement described by the projectors 
n = |0)(0| (a no-click) and III = 1 - II (a click). In 
the rest of this section, we assume that the detector is 
ideal, while ??apd 1 can be taken into account by sub- 
stituting T — ► T]apy)T. 

The probability of no-click of an ideal detector is given 
by P — (0\p\0). It can be determined from the Husimi 
Q-function which provides a phase-space representation 
of the state p: Q(a) is defined as the overlap of p with 
a coherent state \a). The Q-function of a Gaussian state 
with zero mean values of quadratures is a Gaussian func- 
tion centered at the origin, 



Q(r) = 



1 



: exp 



-\r T in + T)-h 



2vrVdet(7 + 7) 
where r = (x,p) and I is the identity matrix. Since vac- 



uum is just a special case of a coherent state, the prob- 
ability of projecting the state (|15fl onto vacuum reads 
P = 47rQ(0) and, on inserting r = in Eq. I)15|l. we 
obtain 



P = 



Vdet( 7 + I) ' 



(16) 



The characterization method works by carrying out 
measurements of the probabilities of no-click Pj for sev- 
eral different transmittances Tj. The covariance matrix 
7j of the state after passing through the beam splitter 
reads 7j = Tj-y + (1 — Tj)I, where 7 is the covariance 
matrix of the input state. On inserting jj into Eq. JT^J, 
we obtain after some simple algebra 



p2= T f det( 7 ) 



■T j (2-T J )Tr( 1 ) + (2-T 1 



(17) 



We thus find that Pj depends on Tj (or, more gener- 
ally, on rjAPuTj) and on the determinant and trace of the 
covariance matrix 7 of the input state. Note that 4/Pj 2 
is a linear function of the two unknown quantities det( 7 ) 
and Tr( 7 ). Thus, measurements of Pj for only two differ- 
ent transmittances simply suffice to determine the trace 
and the determinant, as the system of linear Eqs. i|17|) 
can easily be solved and yields 
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Then, having obtained the determinant and trace of 7, 
we can determine the squeezing properties of the state 
from Eq. © as well as its purity l|l(J|l . 

Dealing with a real world experiment, with unavoid- 
able noises and uncertainties, a more realistic procedure 
would consist in performing the experiment for as many 
transmission values Tj as possible and then trying to 
get the most information from these various measure- 
ments. One possibility to gain information from more 
than two measurements is to implement a maximum- 
likelihood (ML) parameter estimation method (for a re- 
view, see for instance Rcfs. [3 E HI ) This procedure 
provides the values of the parameters Tr(7) and det( 7 ) 
that are the most likely to yield the observed experi- 
mental data. In mathematical terms, this boils down to 
finding the maximum of the joint probability density 



£(Tr(7),det(7)) = n^ Ar "' 



(20) 



which is called the likelihood function of the given cxper- 

(15) 

imental data. Here, Cj denotes the number of photode- 



tector clicks per second for transmittance Tj and N lep is 
the pulse repetition rate. The probability Pj is linked to 
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FIG. 5: Number of photon-detection events per second ver- 
sus average pump power for maximum transmission of the 
variable beamsplitter T = 1. The solid line is a fit following 
Eq.J22l, from which we extracted an estimate of the photon- 
couting detector efficiency 7/apd = 0.84 x 10 -2 ±0.013 x 10~ 2 . 



Tr(7), det(7) and Tj by Eq. (TJJ. Actually, we also have 
to take into account some additional constraints on the 
parameters Tr(7) and det(7). The fact that the covari- 
ance matrix 7 is positive definite and must satisfy the 
generalized Heisenberg uncertainty relation det(7) > 1 
puts the constraints 



1 < det(7) < 



^ Tr( 7 ) 



(21) 



Next section now presents and discusses the results 
of this characterization procedure from its experimental 
implementation. 



V. EXPERIMENTAL RESULTS 

Hereafter, we denote by T only the transmittance of 
the (lossless) variable beamsplitter. Non-unit transmis- 
sions of the spectral and spatial filters and imperfect de- 
tection efficiency are taken into account by an overall 
efficiency parameter t]apd of the APD detection. 

A first step is to estimate this overall efficiency ?7apd 
in order to apply the characterization method. Setting 
the supplementary beamsplitter to a transmittance of 1, 
?7apd can be estimated from the measurement of the 
number of photon-counting events detected per second 
^clicks for different pump powers. In the limit of low 
??apd's, the number of clicks detected per second can be 
approximated as 

Clicks = ^APDiV re p[(# - 1/2) (G + l/G) - 1], (22) 

where the dependence of G and H versus the pump power 
is obtained from the curve fit on the "classical" results 
presented Figs. Eta) and (b), while iV re p = 780.4kHz is 
the repetition rate. With our experimental results (see 



Fig. |SJl and a repetition rate of 780.4 kHz, the fit of A^ c ii c k s 
versus the pump power gives ??apd = 0.84 x 10~ 2 ±0. 013 x 
10~ 2 . This value can be cross-checked with the overall 
efficiency inferred from transmission factors of an intense 
probe beam : the spatial and spectral filters transmit 
respectively 16% and 17% of the probe beam, while the 
APD quantum efficiency is estimated to about 50%, lead- 
ing to an overall detection efficiency of the probe of about 
1.4%. The difference between the latter value and the 
above estimate of ??apd may be explained by slight dif- 
ferences between the modes of the probe (set for maximal 
classical de-amplification) and the squeezed vacuum. 

In our experiment, we used between 4 and 6 different 
settings for the beamsplitter transmittance Tj. For each 
Tj, we performed 100 measurements of the number of 
clicks per second to get a good statistical accuracy on Cj . 
As a result of an appropriate gating of the detection, the 
dark count rate remained reasonably low (about 20s -1 ) 
and was subtracted from the data. 

As shown above, only two different settings of the 
beamsplitter transmittance, T± and T2, are enough to 
extract the value of det(7) and Tr(7) following Eqs. H18|) 
and (|T§)) . Actually, formula ifTHl) indeed leads to an es- 
timate of Tr(7) which is satisfactorily close to the val- 
ues obtained from either homodyne or classical measure- 
ment. However, as far as the determination of det(7) is 
concerned, the formula i|19[l does not give any reliable 
estimate. This results from the fact that, in the exper- 
iment, we have to work with small detection efficiencies 
?/apd < 1 so that small uncertainties on Pi , P% , T%, T2 
have much larger influence on det(7) than on Tr(7). For 
instance, if we take the derivative of Tr(7) and det(7) 
with respect to Pi, we find, 



d Tr( 7 ) 
dPi 

ddet(7) 
rJFi 



?7apdPi 3 ' 
-16 

Iapd 1 



(23) 
(24) 



This shows that in our experimental setup, the determi- 
nant is about 400 times more sensitive to small uncer- 
tainties on Pi than the trace. 

In order to gain information on the determinant of the 
covariance matrix as well as to increase the accuracy of 
the estimate of its trace, we used the full set of mea- 
surements for the different beamsplitter transmittances 
by performing a maximum- likelihood estimation as intro- 
duced in the previous section. The logarithm of the likeli- 
hood function £ given by Eq. (|20J) was computed from the 
measured data, the above estimate of the overall detec- 
tion efficiency t?apd , and the values of the transmittance 
Tj obtained from direct power transmission of the probe 
beam. The global maximum of log(£) was then found by 
brute force numerical search. The experimental results of 
the estimated Tr(7) for several different pump powers are 
shown in Fig.^Ja), and fully coincide with the values in- 
ferred from the classical gain measurements. Out of the 
three trace-estimation procedures, the photon-counting 
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method associated with log-likelihood maximization pro- 
vides the lowest uncertainty on the result. 

Unfortunately, given the low detection efficiency ?7apd 
of our experimental setup, the likelihood function is al- 
most flat as a function of det(7) in the region that is 
allowed by the constraints Q21|) . Consequently, no reli- 
able estimate of the determinant could be obtained from 
our experimental data, the log-likelihood maximization 
method returning essentially a random value between 1 
and (Tr(7)/2) 2 . Similarly, our experimental data only 
provides bounds on the parametric gains G and H given 
the sole knowledge of the trace Tr(7) : 



1 < G < 



1 < H < 



Tr( 7 ) + y / TV(7T^4 
Tr( 7 ) - VTr( 7 ) 2 - 4 
Tr( 7 ) + 2 



1/2 



(25) 
(26) 



Thus, the dashed lines in Figs. E{a), (b), (c), and^Jb) 
take into account the fact that no estimate of the deter- 
minant better than the bounds l|21|l could be obtained 
by the photon-counting method given the estimate of the 
trace. We also tried various other numerical methods - 
such as least squares inversion - but neither provided a 
reliable estimate of det(7). 

Some better insight on the intrinsic difficulty to get an 
estimate of det(7) can be obtained by rewriting Eq. l(T7|) 
as 



= (det(7)-Tr( 7 ) + l)7 ? i PD T; 
+ 2(Tr( 7 )-2)ry APD T j +4. 



(27) 



It becomes clear that the determinant is linked to the 
second-order dependence of PJ 2 in the transmittance, 
while the trace can be directly obtained from the lin- 
ear dependence of PJ . The basic difficulty to estimate 
det(7) results from the fact that the relevant informa- 
tion is hidden in terms of order (tjapdT) 2 , which are very 
small for our experimental data given the low values of 

?7APD ■ 

One could then try to increase the overall APD de- 
tection efficiency ?7apd by releasing cither the spatial or 
spectral filtering conditions. However, from an exper- 
imental point of view, this does not seem realistic for 
several reasons. 

First, we would move to a region where the physics be- 
come multimode, which is clearly outside the framework 
of the developed model. In principle, the photon count- 
ing method allows one to check whether the single-mode 
description of the experiment is appropriate or not. If 
only a single mode is detected, then P~ 2 should be a 
quadratic polynomial in i]app,T, cfr. Eq. I|27l) . More 
generally, if the detector effectively registers light from 
N modes in a Gaussian state, then P~ 2 becomes a poly- 
nomial of 2iV-th order in t/apbT fl3| . So, after measuring 
P as a function of t/apdT one could perform a fitting to 



determine the minimal number of modes N that is neces- 
sary for the description of the observed signal. However, 
a successful application of this technique would require a 
very high precision in the measurement of P and a high 

??APD- 

A second problem with removing the spatial and/or 
spectral filters is that we would loose any possibility to 
cross-check our results with classical parametric gain or 
homodyne measurements. Last, even in the case of no 
spatial filter and 10 nm spectral filter, the overall APD 
detection efficiency will remain low given our experimen- 
tal setup, and we do not expect to gain much according 
to our numerical simulations of the constraints on the 
global efficiency presented below. 



VI. NUMERICAL SIMULATIONS 

We have seen that the low APD detection efficiency 
t?apd precludes a reliable estimate of det(7) via the 
photon-counting method. It is thus important to deter- 
mine the efficiency t?apd that should be attained in order 
to be able to estimate det(7) with acceptably small er- 
rors. More generally, it is interesting to investigate the 
dependence of the estimation errors on ?7apd- For this 
purpose, we have carried out extensive numerical sim- 
ulations of the experiment for several values of ?7apDj 
the other parameters of the simulation being chosen in 
accordance with the experimental values. In particu- 
lar, we have assumed a measurement repetition rate 
N Iep = 780.4 kHz and a total measurement time t = 100 s 
for each transmittance Tj . The total number of measure- 
ments for each Tj is then given by iVtot = N Iop t. We 
have further assumed that measurements were carried 
out for four different transmittances T\ = 1, T% = 0.75, 
T 3 = 0.5 and T4 = 0.25, and we used the experimentally 
obtained values det(7) = 1.156 and Tr(7) = 2.321 as a 
typical example (corresponding to a pump average power 
ofl.21mW). 

The determinant and the trace of 7 were estimated 
from the simulated experimental data with the help of 
the maximum-likelihood technique described in the pre- 
ceding section. Since the ML estimator is generally bi- 
ased, we define the deviation of the estimate from the 
true value as 



^dct = ((det( 7 ) cst 



det(7)truo) 2 ), 



"it 



= <(Tr( 7 ) cst - Tr( 7 ) true ) 2 ), 



(28) 



where () indicates averaging over an ensemble of exper- 
iments. In practice, we simulated 1000 times the whole 
experiment, from data acquisition to the ML estimation, 
and we then calculated (|28|l by averaging over the en- 
semble. Since the total number of measurements iVtot 
was very large, we approximated the binomial distribu- 
tion of Cj by a normal distribution with the same mean 
and variance. 

Besides the statistical fluctuations of Cj and the in- 
trinsic difficulty of estimating det(7) at low detection effi- 
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FIG. 6: The dependence of the variance <7d e t °f the estimation 
of det(7) on the detector efficiency t^apd in case when the 
parameters Tj and t^apd are known exactly (circles) and when 
the experimental uncertainties of Tj is 0.5% and the relative 
uncertainty of ?7apd is 1% (squares). See text for further 
details. 
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FIG. 7: Results of the estimation of det(7) from simulated 
photon counting measurement assuming overall detection ef- 
ficiency tjapd = 50%. This figure has been scaled so as to 
be easily compared to the experimental results presented on 
FigHb). 



ciencies, other factors contribute to the estimation errors, 
namely the uncertainty in the knowledge of Tj and ?7apd ■ 
To isolate the errors stemming from low ?7apd, we have 
first assumed that all parameters Tj and t?apd are known 
precisely, hence the statistical fluctuations of Cj are the 
only source of errors. The resulting a^ct is plotted as cir- 
cles in Fig. For very low t^apd, the estimates of det(7) 
are randomly distributed in the interval [1, (Tr(7)) 2 /4] 
and crdet ~ |(Tr 2 (7)/4— 1). The estimation error rapidly 
decreases as ?7apd grows, and our numerical simulations 
reveal that a reliable estimate of det(7) with ovjet < 10~ 2 
could be obtained for ?7apd > 15%. 

The uncertainties of Tj and ?7apd significantly increase 
the estimation error for higher ?7apd • We have performed 
numerical simulations taking into account that TVs are 
known with an uncertainty of 0.5%, and the relative un- 
certainty of 7?apd is 1% which corresponds to the actual 
experimental situation. The resulting cr^ct is plotted as 
squares in Fig. El We observe that Cdot is much higher 
than in the previous case, except for the region of very 
small ?7apd- To obtain a satisfactorily accurate estimate 
of det(7) with a w 2 x 10~ 2 , we need t^apd > 50%. 

In order to demonstrate that ?7apd = 50% is indeed 
sufficient for the whole range of values of the pump power, 
we have simulated the results of an experiment at t^apd = 
50% for the same values of the pump power as in Figs.|3 
and 01 The results are given in Fig. [3 which shows the 
mean estimated values of det(7) as well as the resulting 
error bars. We find that these estimates are in very good 
agreement with the true values used in the simulation. 

Finally, note that our numerical simulations also con- 
firm that the estimate of Tr(7) is very accurate: we have 
found that orr < 10~ 2 even for t^apd as low as 1%. 



VII. CONCLUSIONS 

In this paper, we have shown that direct photon- 
counting detection can be used, instead of homodyn- 
ing, to evaluate the squeezing and purity of an arbitrary 
single-mode gaussian state. For the rather generic states 
that we considered, the trace of the covariance matrix can 
be accurately determined, even with an overall detection 
efficiency r\ in the percent range, while its determinant 
(related to the state purity) requires a much higher 77, 
typically around 50%. 

In principle, such efficiencies are well within the reach 
of silicon photon-counting avalanche photodiodes, but an 
important problem remains : most sources do not emit 
single- mode gaussian light, but rather multimode light. 
This is not a problem when a homodyne detection is used, 
because the local oscillator acts as a very efficient single- 
mode filter. On the other hand, a photon counter de- 
tects photons in any mode. Therefore, detecting a good 
approximation of a single mode state requires appropri- 
ate spatial and spectral filters, respectively obtained from 
pinholes and diffraction gratings. Unless a special effort 
is made, these filters will have a low overall transmission 
(a few percent in our experiment), and thus the direct 
detection method will fail to determine the state purity. 

In principle, there are various ways for improvement, 
which are open for further experimental work. Ideally, 
the source itself should emit single mode light, which 
might be obtained by appropriate phase-matching condi- 
tions in a x^ 2 ) non-linear crystal. On the filtering side, in- 
tcrfcromctric multidielectric filters provide transmission 
values which are much higher than standard slits and 
grating set-ups. A combination of these various tech- 
niques will be probably needed to reach the high overall 
efficiencies needed for many potential applications. 

As a conclusion, it appears that a broad variety of 
techniques is available to characterize quantum contin- 
uous variables, and that these methods will certainly 
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continue to develop for applications in quantum cryp- 
tography, quantum communications, and possibly quan- 
tum computing. Perhaps the most appealing application 
of the photon-counting method is the direct determina- 
tion of the entanglement of two-mode Gaussian states by 
measuring only the purity of the two-mode state and the 
marginal purities of the single-mode states on each side 
[l3L |2]| . All these purities can be determined with the 
photon counting method using only local measurements. 
The distinct feature of this approach is that no interfero- 
metric stability is required if one is dealing with squeezed 
vacuum states, which is the case in many experiments. 
This may be an important advantage in the character- 
ization of entanglement distribution over long-distance 



continuous- variable quantum communication networks. 
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